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1 Foamlike Patterns in Cosmic Structure 

Probing cosmic large scale structure on the basis of X-ray observations puts par- 
ticular emphasis on the densest regions within the global matter distribution, the 
rich clusters of galaxies. Understanding the relationship between the cluster dis- 
tribution and the underlying matter distribution is therefore a key element in any 
assessment of cosmic structure on the basis of samples of X-ray selected clusters. 
By now, the foamlike arrangement of matter and galaxies on Megaparsec scales 
has become a well-established feature of the cosmic matter distribution, consist- 
ing of an assembly of anisotropic elements, filaments and walls of various sizes, 
surrounding large underdense void regions and sprinkled with wholly or partially 
virialized dense clumps of matter, varying in size from rich clusters of galaxies 
down to small groups of a few galaxies. It is one of successes of gravitational 
instability cosmic structure formation theories to find that this pattern of walls, 
filaments and voids appears to be the generic outcome of these scenarios. 

A major obstacle in quantifying this quintessential aspect of cosmic structure 
is the lack of a systematic insight into the dynamical and statistical aspects of cel- 
lular geometries as well as the absence of a readily available and well-established 
mathematical machinery to evaluate and compare observations and simulations. 
Stochastic geometry - the branch of mathematics concerned with nontrivial geo- 
metrical concepts involving stochastic behaviour of one or more of their character- 
istics - may be expected to contribute significantly to further such understanding. 

2 Voronoi Tessellations 

The canonical example of a stochastic geometrical model for a cellular division 
of space is that of Voronoi tessellations. This space filling network of convex 
polyhedra offers a surprisingly realistic and versatile representation of the char- 
acteristics and features of the foamlike or cellular spatial arrangement of matter 
in the Universe. In short, it is defined through a spatial distribution of nuclei 
Each nucleus corresponds to one Voronoi cell, which comprises that part of space 
closer to this nucleus than to any of the other. The walls and edges forming 
the polyhedrons' surfaces are identified with the wall-like and filamentary su- 
perclusters in the galaxy distribution, the vertices with and massive clusters of 
galaxies, while the interior of the Voronoi cells corresponds to the large void re- 
gions barren of galaxies. The morphology of artificial galaxy distributions set up 
within this network bears a striking resemblance to that observed in the more 
complicated circumstances of the real Cosmos or that of the artificial reality of 
computer simulations of structure formation. Figure 1 shows a conglomerate of 



Figure 1: A clump of Voronoi polyhedra from a tessellation of 1000 cells. 
One cell is surface-shaded, the surrounding neighbouring cells are depicted 
in wire-frame fashion. The corresponding vertices are indicated by solid 
dots. Distribution of all Voronoi vertices within a cubic volume. Courtesy: 
Jacco Dankers. 

several neighbouring Voronoi cells, with one of the cells by shaded surface and a 
wire-framed representation in the case of its neighbouring cells. The solid lines 
are the edges (filaments) in the network, the dots (coloured red) are the vertices 
of the Voronoi tessellation. In principle a Voronoi tessellation could be regarded 
as a mere geometric toy model, a heuristic geometrical description of a substan- 
tially more complicated reality. However, a detailed assessment of the role of 
voids in structure formation provides a physical ground for them representing re- 
ality, whether the observed or simulated one, in a more subtle way [HQ. This 
is certainly reinforced by their striking resemblance to the observed or simulated 
foamlike appearance of the galaxy and matter distribution. 

3 Voronoi Vertex Clustering 

Even in the limited setting of Figure 1 it is evident that the vertex distribution is 
not a random Poisson distribution. The full spatial distribution of Voronoi vertices 
in the full cubic volume (righthans frame fig. 1) clearly involves a substantial 
degree of clustering. This impression of strong clustering, on scales smaller than 
or of the order of the cellsize Ac, is most evidently confirmed by their two-point 
correlation function £(r). Not only can we discern a clear positive signal but - 
surprising at the time of its finding on the basis of similar computer experiments || 
- out to a distance of at least r ~ 1/4 Ac the correlation function appears to be 
an almost perfect power-law, 



£(r) = (r /r) 



(1) 



with a slope 7 ~ 1.9 — 2.0. Its amplitude, traditionally expressed in terms of the 
"clustering length" r a , at which £(r ) = 1, has a value r » 0.29 Ac- Beyond 
this range, the power-law behaviour breaks down and following a gradual decline 
the correlation function rapidly falls of to a zero value once distances are of the 
order of the cellsize. However, rather than a characteristic geometric scale, r is 
more a measure for the "compactness" of the spatial clustering, set mainly by the 
small-scale clustering. A more significant scale within the context of the geometry 
of the spatial patterns in the density distribution is the "correlation length" r a , 
the scale at which f (r a ) = . As a genuine scale of coherence, it is more relevant 
to the morphology of the nontrivial spatial structures we seek to study. Beyond 
r a the distribution of Voronoi vertices is practically uniform. 

If we interpret the clusteringlength r ~ 20/i _1 Mpc, usually found for sam- 
ples of rich clusters of galaxies H, within the context of a Voronoi tessellation 
it would imply a cellsize of Ac ~ 70/i _1 Mpc. Although the two-point cluster- 
cluster correlation function reproduced by the Voronoi vertices fits very well to 
the function yielded from the observations, the large cell size A c may be a com- 
plication. It is surely well in excess of the 25/i -1 -35h -1 Mpc size of the voids in 
the galaxy distribution. Moreover, also within the Voronoi concept itself it would 
conflict with the clustering of objects dwelling in the walls and filaments of the 
same tessellation framework. Clustering analysis of such configurations || reveals 
that the two-point correlation function of galaxies confined to the walls - as well 
as for those confined to the edges - also displays distinct power-law behaviour 
at sub-cellular scales. The involved clustering length, however, is different from 
that of the vertices in the same framework. For the wall galaxies it is but half 
the value of that of the vertices, r m>Q ~ 0.14 Ac||,[jl0]]. If 7" Wi0 is identified with 
the galaxy-galaxy clustering length of r gj0 ~ 5h~ Mpc, this would yield a cellsize 
of ~ 35/i -1 Mpc. The latter is suggestively similar to the size of the actually ob- 
served voids, which may be a tantalizing hint for a profound relationship between 
the clustering length r and the typical cellular scale of the cosmic foam network. 

Adressing this apparent inconsistency between vertex and wall clustering we 
first observe that the vertex correlation function of eqn. (1) concerns the whole 
sample of vertices, irrespective of any possible selection effects based on one or 
more relevant physical aspects. In reality, it will be almost inevitable to invoke 
some sort of biasing through the defining criteria of the catalogue of clusters. In- 
terpreting the Voronoi model in its quality of asymptotic approximation to the 
galaxy distribution, its vertices will automatically comprise a range of "masses" . 
Neglecting the details of the temporal evolution, we may assign each Voronoi ver- 
tex a "mass" equal to the total amount of matter that ultimately will flow towards 
that vertex. Applying the "Voronoi streaming model" || as a reasonable descrip- 
tion of the clustering process, it is reasonably straightforward if cumbersome to 
calculate the "mass" or "richness" A^v of each Voronoi vertex by pure geomet- 
ric means[jn^]. It concerns the volume of a non-convex polyhedron centered on 
the Voronoi vertex, with related Voronoi nuclei as one of the polyhedral vertices. 
These nuclei are the ones that supply the Voronoi vertex with inflowing matter. 
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Figure 2: Voronoi tessellation scaling. Left insert: Two-point correlation 
function £(r) Voronoi vertices. The clustering length r Q (left) and the ratio 
r&/r (right) as a function of average vertex distance. 



Evidently, the "mass" is larger for vertices on the surface of large Voronoi cells. 

The vertex samples in our study consist of all Voronoi vertices with a richness 
equal to higher than a limiting value, the "sample richness" -Ms[0- Subsequently, 
we determined the correlation function characteristics for each of the vertex sam- 
ples. Assessing their behaviour as a function of the average distance Av between 
the sample vertices, for samples encompassing from 10% up to 100% of the total 
number of Voronoi vertices and thus for Av ~ 0.5 — 1.5 x Ac, revealed a remarkable 
and tantalizing scaling of the clustering phenomenon [^o| . 

All subsamples of Voronoi vertices have a two-point correlation function that 
out to a certain range retains a power-law behaviour almost exactly similar the 
one of the full sample of vertcies. No significant difference in power law slope 
between the various subsamples can be discerned, all have 7 ~ 1.8 — 2.0. On 
the other hand, both clustering length r and correlation length r a do display a 
systematic dependence on sample richness. Both r and r a increase proportionally 
to the sample richness, and hence the average vertex distance Av in the samples ! 
The more massive the vertices are, the more strongly they cluster. In other words, 
in the line of the random field "peak biasing" scheme pi, we here find a purely 
geometrically based biasing scheme. Interesting in this respect is to observe that 
the increase of rv,o is perfectly linearly proportional to the mutual vertex distance 
Av in the sample (Fig. 2, left frame). This suggests that Voronoi vertex clustering 
is a perfect realization of the clustering scaling description proposed by Szalay & 
Schramm Q. It also adheres to the increasing level of clustering that selections 
of more massive clusters appear to display in large-scope N-body simulations || 
although there are telling differences in detailed behaviour. 

For our purpose, even more significant is the uncovered scaling of the cor- 



relation length rv,a, similar in character to that of ry,o- The increase of r a also 
turns out to be almost exactly linearly proportional to the average vertex distance. 
The repercussions are manyfold. It means that the selected vertices still have a 
strong positive correlation at scales where the poorer samples do not possess any 
clustering. Moreover, samples with more massive clusters are expected to have a 
clustering that extends out further than that of their poorer brethren. This may 
be a significant observation in the light of the finding that samples of rich galaxy 
clusters seem to have a positive correlation at scales of tens of Megaparsec, quite 
in excess of scales with appreciable galaxy clustering. Finally, the implied con- 
stant ratio between clustering and correlation length, r a /r » 1.86 (Fig. 2, right 
frame), implies a perfect self-similar scaling of the Voronoi vertex distribution. 
The complete correlation function £ s (r) of each selected subsample s, and not 
just the part in the power-law range, is a self-similar mapping of an elementary 
function £ e i scaled by means of a characteristic lengthscale parameter L 3 . 

4 Bias and Cosmic Geometry: Conclusions 

The above results form a tantalizing indication for the existence of self-similar 
clustering behaviour in spatial patterns with a cellular or foamlike morphology. It 
might hint at an intriguing and intimate relationship between the cosmic foamlike 
geometry and a variety of aspects of the spatial distribution of galaxies and clus- 
ters. One important implication is that with clusters residing at a subset of nodes 
in the cosmic cellular framework, a configuration certainly reminiscent of the ob- 
served reality, it would explain why the level of clustering of clusters of galaxies 
becomes stronger as it concerns samples of more massive clusters. In addition, 
it would succesfully reproduce positive clustering of clusters over scales substan- 
tially exceeding the characteristic scale of voids and other elements of the cosmic 
foam. At these Megaparsec scales there is a close kinship between the measured 
galaxy-galaxy two-point correlation function and the foamlike morphology of the 
galaxy distribution. In other words, the cosmic geometry apparently implies a 'ge- 
ometrical biasing" effect, qualitatively different from the more conventional "peak 
biasing" picture tefl. 
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